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GRB 080916C AND GRB 090510: THE HIGH-ENERGY EMISSION AND THE AFTERGLOW 
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ABSTRACT 

We constrain the physical composition of the outflows of GRBs 080916C and 090510 with the prompt 
emission data and find that the former is hkely magnetic, while the latter may be baryonic. The X- 
ray and optical afterglow emission of both GRBs can be reasonably fitted using the standard external 
shock model but the density profiles of the circum-burst medium are different. We also propose a 
simple method to estimate the number of seed photons supposing the GeV afterglow photons are due 
to the inverse Gompton radiation of external forward shock electrons. The seed photons needed in 
the modeling are too many to be realistic for both events. The synchrotron radiation of the forward 
shock seems able to account for the GeV afterglow data. 

Subject headings: Gamma rays: bursts — Radiation mechanisms: non-thermal — X-rays: general 



1. INTRODUCTION 

Gamma-ray bursts (GRBs) are the most luminous 
explosions in the universe. They feature extremely rela- 
tivistic outflows with bulk Lorentz factors ^ 10^"'^ and 
isotropic energies of lO***"^^ erg. Though their cosmo- 
logical origin as well as the relativistic movement have 
been flrmly established, the radiation mechanism and 
the outflow composition are still uncertain (Piran 199^: 
iZhang fc^ cszaros 2004). It is widely believed that the 
high-energy emission of GRBs can shed light on these two 
fundamental issues (see Fan & Piran 2008 for a review). 
For example, a distinct GeV-TeV spectrum excess can 
be taken as an indication evidence of a baryonic outflow 
and a radiation process in addition to synchrotron (e.g., 
inverse Compton scattering) will be needed, while the 
absence of such a component in most spectra may favor 
the magnetic outflow model. Recently, the Fermi col- 
laboration has released their observation d ata of GRBs 
080916C and 090510 (|Abdo et al. I I2009allbh . In this 
work, we examine the origins of these prompt and after- 
glow GeV emission. The work is structured as follows. 
In Section 2, we discuss the origin of the prompt GeV 
emission and the corresponding constraint on the phys- 
ical composition. In Section 3 we employ the standard 
external forward shock model to interpret the X-ray and 
optical afterglow data. In Section 4, we investigate the 
origin of the afterglow GeV emission. Our results are 
summarized in Section 5 with some discussion. 



2. PROMPT GEV EMISSION OF GRBS 080916C AND 090510 



GRB 080916C was a lone 



^as a long 
Tgo ~ 66 s (lAbdo et ani2009ali , - 
z ~ 4.35 ± 0.15 (IGreiner et al.ll2009bfl 



burst with a duration 
and was at a redshift 
A few hundred 



high-energy photons have been detected by the large area 
telescope (LAT) onboard the Fermi satellite and three 
of them are above 10 GeV. The joint analysis of the 

Electronic address: |yizhong@nbi.dk (YZF)| 



LAT and Gamma-ray Bursts Monitor (GBM) data sug- 
gests a featureless Ban d spectrum in the energy range 
8 keV - 10 GeV (Abd o et al. 1 l2009a[ ) . A straightfor- 
ward interpretation of the spectrum is the synchrotron 
radiation of internal shock electrons. Such an interpreta- 
tion, if correct, demands a very large bulk L orentz factor 
Fj ~ 1 0'^ of the emi tting/shocked region (|Abdo et al. I 
I2009at iGreiner et all [20"09b) . In the internal shock sce- 
nario, the fast shells should move faster and the corre- 
sponding bulk Lorentz factor should be Ff ~ 5Fi other- 
wise the internal shock efficiency will be too low to match 
the observations (e.g.. iPiraniil999 ). The photosphere ra- 
dius of the fast shells is i?ph ~ 5 x 10^ cm L54F^3y 
(jPaczvhski II1990I ). where L is the total luminosity of the 
outflow^. On the other hand, for a baryonic shell we 
have Ff < i?ph/i?o 5 x IQ^ Lr^iTJ^^R^l ()Piranlll999f l. 
where i?o > 10^ cm is the size of the central engine. So 
the shell becomes transparent at the late stage of its ac- 
celeration. As a result, the thermal radiation of these 
shells will be too strong to be effectively outshone by 
the internal shock non-thermal emission, in disagreement 
with the data (Fan 2009; see Zhang & Pe'er 2009 for the 
other approach). Hence we would not discuss the stan- 
dard/unmagnetized internal shock model for this burst. 

An interesting possibility is that the prompt emission 
has a very soft MeV-GeV spectrum and the GeV pho- 
tons are due to the synchrotron radiation of the external 
forward shock (Kumar & Barniol Duran 2009). Here we 
outline a few potential challenges of such a model. (1) 
In the forward shock model, the variability of the radia- 
tion is determined by the angular timescale Tang, which 
is '--^ t as long as the edge of the emitting region is invis- 
ible (Piran 1999). So the light curve should be smooth. 
The variability shown in the LAT data then disfavors the 
forward shock emission model. (2) For the initial outflow 

^ In this work we adopt the convenience Qx = in units 

of cgs. 
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expanding into the wind medium (see Section 3 for the 
medium identification), strong reverse shock may form. 
The bulk Lorentz fact or of the s hocked medium will be 
almost a constant f C hevalier fc Li 2000). A strong re- 
verse shock exists till t ~ 190/2. In such a phase, we 
have the magnetic field strength B oc t~^, the maxi- 
mum specific flux -Fl/^max oc t'^ , the typical synchrotron 
frequency t'm oc and the cooling frequency t. 
Hence the synchrotron radiation flux in LAT band can 
be estimated as -Flat oc Fl^.max^'m^^^^^i^y^ oc t(2-p)/2 
for hvc < 100 MeV, inconsistent with the observation. 
Where p is the power-law distribution inde x of the ac- 
cele rated electrons at the shock front (see iXue et alTI 
12009;. for extensive discussion). Since the reverse shock 
emission has not been detected in most GRBs and it is 
not clear whether the model suffers some disadvantages, 
we do not take the current temporal inconsistence as a 
conclusive argument. (3) To reproduce the prompt spec- 
trum, the forward shock emission at t ~ 10 s should have 
hi^m > 300keV. At such early time, the synchrotron self- 
Compton radiation is in extreme Klein-Nishina regime 
and the Compton parameter Y ^ 0. With proper pa- 
rameters, Uc can be comparable to z/m- So the sub-MeV 
spectrum can be F^, oc i/^/^, steep enough to be consis- 
tent with the data. However, if ~ 10'^°{t/10)~^^^ Hz, 
the XRT light curve will be F^^ cx for i < 10^ s and the 
optical light curve will be F^^_^^ oc t° for t < 10^ s. These 
behaviors are very unusual and have not been detected 
in other GRB afterglows so far. The lack of observation 
of early afterglow of GRB 080916C, however, hampers 
us to test the model. 

If the prompt high-energy emission of GRB 080916C 
was from the soft gamma-ray emitting region, a plausible 
origin of the GeV photons is the synchrotron radiation of 
electrons accelerated in magnetic energy dissipation of a 
Poynting-flux dominated outflow (jZhang fc Pe'erj[2009D . 
A disadvantage of such a scenario is the difficulty of re- 
producing the hard low energy spectrum (jFan II2009I ). 

G RB 090510 was a short burst at a redshift z ^ 
0.903 (jAbdo et al. 11200 9b'). The high-energy emission 
is much more intense than that of GRB 080916C and 
shows some variability, which disfavors the external for- 
ward shock model. In the time interval 0.5 — 0.6 s, 
the sub-MeV spectrum is very hard but the h igh en- 
ergy spectrum is very soft (Ab do et al. Il2009bf ). possi- 
bly dominated by the photosphere emission of the bary- 
onic shell^. In the time interval 0.5 — 0.8 s, the high en- 
ergy spectrum gets harder and harder but the "thermal" - 
like MeV component is still evident. GeV emission is 
naturally produced in the IC scattering of the "photo- 
sphere" photons by the shocked electrons. The pho- 
tosphere radius is ~ 6 x 10^^ cm L54T~^^, where Tgh 
is the bulk Lorentz factor of the shell. The inter- 
nal shocks take place at a rather larger radius i?^ ^ 
T?cSt/{l + z) 1.5 X 10^5 cm Tl^{St/0.1 s), where dt 
is the detected variability timescale of the prompt emis- 
sion. In the comoving frame of the emitting region the 

The temperature of the initial shell is Tobs ~ 10 MeV [(1 + 
2)/2]~^L5^^i?Q , matching the data if Rq ~ 10^ cm. Consider- 
ing the un-magnetization nature of the outflow, such a small Ro 
indicates a black hole as the central engine. The outflow was likely 
launched via the neutrino-antineutrino annihilation process. 



seed/photosphere photons are moving along the radial 
direction and are highly anisotropic. In such a case, the 
strongest IC radiatio n is from an angle ~ l/Fi relative 
to the line of sight (|Fan fc PiranI [2006b) . The arrival 
of the GeV photons will be delayed by a time ~ dt 
and the GeV radiation duration will be extended, in 
agreement with the observation. Below we show how 
to reproduce the high energy spectrum F,^ oc ;y~0-54 jj^ 
time interval 0.8 — 0.9 s. If the cooling of the electrons 
is dominated by the prompt soft gamma-rays with a 
luminosity L^, the cooling Lorentz factor can be esti- 
mated by 7c^ic ^ 5 £,753 gi?-. i.^rf g ilFan fc PiranI 120081 ). 
Here we do not take ^ 10^^ erg s~^, the lumi- 
nosity of the simultaneous soft gamma-ray emission, 
since in the photosphere-internal shock model the ar- 
rival of the upscattered photons is delayed, as already 
mentioned. The corresponding IC radiation frequency 
ec,ic ~ Tc'ic^^P 25 MeV (E^/l MeY)L-% ,Rl,,,:,rl,, 
where Ep is the typical energy of the seed photons. 
On the other hand 7,n i « Cc i{mp/me)(Tsh — l)/3 ~ 
100(ee,i/0.5)[(Fsh - l)/d.3] for p ~ 2.5, where F^h is the 
parameter denoting the strength of the shocks. The cor- 
responding IC radiation frequency is Em.ic ~ 7m i-^p ~ 
10 GeV (7,„,i/100)2(£'p/l MeV). The spectrum in the 
energy range ~ 10 McV — 10 GeV is F^^ oc v^^^^, con- 
sistent with the data. We note that in the time interval 
0.9 — 2 s, the soft gamma-ray emission is very weak while 
the GeV emission is still strong. These delayed GeV pho- 
tons may be produced by the IC scattering of the soft 
gamma-rays by the electrons accelerated by the reverse 
shock or by the shocks generated in the collision of the 
late time {t > 0.5 s) outflow with the precursor outflow. 

3. THE AFTERGLOW OF GRBS 080916C AND 090510 

GRB 0809 16C. Swift XRT started to observe this 
source at about 17 hr after the Fermi trigger. In our 
data analysis, the X-ray light curve can be fitted by a 
single power-law F^^ oc i-i-30±0.07 g ^ x 10^ < t < 
1.3 X 10^ s and the XRT spectrum is F^ oc :y-0-50±o.i6_ 
The earliest optical/infrared observation started aX t ^ 
26.7 hr after the burst. The optical/NIR light curve can 
be well described by F^^^^ oc t-i-40±o.05^ ^^Yie optical 
to X-ray spe ctrum is consistent w ith a single power law 
F^ oc v-°-^^ ()Greiner et al.ll2009bf) . These facts suggest 
that the optical to X-ray afterglow emission is within 
the same regime. In the standard external shock model 
(e.g., Zhang fc Meszaros 2004), the slow cooling spectrum 
takes the form F^, oc and the decline should be 

either t3(i-p)/4 (ISM) or t(i-3p)/4 (wind medium). One 
can see that the X-ray and optical afterglow data are in 
agreement with the wind medium model for p ~ 2.2 (see 
also Zou et al. 2009). 

Assuming a GRB efficiency e ~ 0.2, the isotropic- 
equivalent kinetic energy of the outflow is i?k ~ 4 x 
gj.gg ^jjg wind case, the eqt iations that gove rn 
the forward shock emission are (e.g.. lYost et~all[2003l) 

v^, « 1.3 X 10^4 Hz el_,e'J,^ClEl% ,{l + zfl%^\ 

(1) 

« 1.7x10" Hz .//X%'5 6^:'(l+-)-'/'^14^l+n'^ 

(2) 
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F.,„,ax « 100 mJy + ^f'^^Dl^.,. 

(3) 

where Cp = 13{p - 2)/[3{p - 1)] for p > 2.05, 
Ah, = (M/IO-^Mq yr-i)[ww/(10^ cm s"!)]-! is the 
wind parameter Vy, is the speed of t he wind, M is the 
mass loss rate (jChevalier fc Lil [2000l ) . and Y = [-1 + 

v/l + 4?7?7j,„ec/eB]/2, ~ mm{l, and 
77j^„ is the factor reflecting the importance of the Klein- 
Nishina correction (see the Appendix A of iFan &: PiranI 
(|2006aD for the expression). 

Since v^n decreases with time while Vc increases with 
time, the current afterglow data suggest that Vyn{t = 
10^ s) < J/opt/iR and z^c(i = 6 X 10* s) > z/^ 10^** Hz, 
i.e.. 



,2 1/2 



<5, 



.-3/2.-2 



(1 + > 7 X 10^ (4) 



At t ~ s, th e Ks band flux is ~ 3 x 10"^ Jy 
(|Greiner et al.|[2009bf ) , which gives us another constraint 



Lie^M, - 7 X 10" 



(5) 



Substituting Y ^ eo/50eB (due to the slow cooling 
and the Klein-Nishina correction) in Equations (4) and 
(5), we have A„ > IQ-^e^ cb > lO-^e;;^^ Though 
the shock parameters cannot be uniquely determined, 
we see that the "reasonable" parameters (ee, cbi ^*) ~ 
(0.1, 2.5 X 10~^, 0.01) can reproduce the data. 

GRB 090510. In our data analysis, before and 
after the break at th ~ 676(1 + z) s the X-ray 
declines are t-0-72±o.08 ^-i.89±o.06^ respectively. 

The X-ray spectrum can be reasonably fitted by 

oc i/-0-63±o.06_ reduced the UVOT data in 

a standard way with t he aid of reduction threads at 
Ihttp:// www . swift . ac . uk / U VO T . shtml , The combined V- 
band and white light curves sliow a rise since the begin- 
ning of UVOT observation to a peak around 1000 s after 
the BAT trigger, which is followed by an apparent decay 
leading to the optical flux lower than the threshold of 
UVOT quic kly. Our results ar e generally in agreement 
with that of iDe Pasquale et all |2009) . Within the stan- 
dard external shock model, the above data are roughly 
consistent with a slow cooling ejecta expanding into the 
ISM for p ^ 2 while the break can be interpre ted as 
the jet effect (|Piranlll999l : lZhang fc Meszarosll2004f ). The 
slowly rising optical emission may suggest that the ob- 
server's frequency is below v^. In the ISM case, the 
equations that govern the forward shock emission are 
(e.g., Sari et al. 1998; Yost et al. 2003) 

z/, « 5.2 X IQi* Hz i^k^sf es'^Wll + + Y 



z... = 7.0 X 101^HzE^/5Vb(^4<-iC^(1 + z)'/'t3.^/^(7) 



F.,niax = 2.7 X 10-3 Jy (1 + 7)\ill^^^^e]^^_^^i,^^inl'\^) 

please note that we have Cp ~ 0.23 for p ^ 2. 

The conditions that Vc{t — 1284 s) > v^^ v^{t ~ 
1000 s) ~ 5 X 10^4 Hz and i^^.^ax > 1 x lO"* Jy 
(|De Pasquale et al.ll2009[ l yield 



^i:5f^B'-W(i+n"'>o.2, 



El%4U<-i - e'J'-,EkMnl^' > 0-02. (10) 

The parameters (£'k,54, eB,-4, Ce.-i, "■o) ^ 
(1, 1, 7, 0.01) satisfy the above constraints (note that 
Y <C \/ Ec/eB thanks to the Klein-Nishina correction). 
The jet break time — 1284 sec suggests a half- 
opening angle 0j = 6 x W^^t^^,iE^l4^f^^Jo,Ynl^l2- So the 
true gamma-ray energy released is fi'^jct — S?E^/2 — 
2 X lO^'^ergs, where ~ 1.4 x 10^"^ erg is the isotropic- 
equivalent gamma-ray energy. 

4. THE HIGH ENERGY AFTERGLOW EMISSION 

4.1. IC scattering in the forward shock region? 

If the high energy afterglow is due to the IC radi- 
ation of the forward shock electrons, there is a simple 
method to estimate the number of seed photons, regard- 
less of their origin (either the late prompt emission from 
the central engine or the synchrotr on radiat ion of the for- 
ward shock electrons). Following IFan fcT'iranI l|2006bf) . 
the possibility of one seed photon being scattered (i.e., 
the optical depth) in the forward shock region can be 
estimated as 

risM - 4.2 X 10-« El%nl/*tl/'[{l + z)/2]-^/\ (11) 

Twind - 7.3 X 10-6 ATE;li%'/'[{l + z)/2]i/2, (12) 

respectively. With the parameters derived for GRBs 
080916C and 090510, we have 

rwind(080916C) ~ 10-^ Al%E-lilt-'/\ 

tism(090510) ^ 7 x 10-^° Eli^^n'J^tl^^ 
respectively. 

If the detected high energy afterglow photons are 
indeed the IC radiation of the forward shock electrons, 
the number flux of the seed photons will be 



seed 



>100MeV 



/r. 



(13) 



For 



GRB 080916C, in the time inter- 
val - 100 - 1400 s (i.e.. At = 1300 s), 
-F>iooMcV 7 X 10-6 ph cm-2 s'^ (Abdo et al. 2009a), 
so the number of total seed photons is 



(TT^ 



Ati^seed ~ 6 X 10' 



64 



(14) 



(9) 



If most seed photons are in the X-ray band, the total 
nergy will be ~ 10^^ gj-g, which is too large to be real- 
'stic. If the seed photons are mainly in optical/infrared 
band, the total energy will be ^ 10^^ erg. Though bright 
infrared/optical flare can be produced in the afterglow 
phase by the prolonged activity of the central engine 
(for example, t he infr a red flare det ected in GRB 080129; 
iGreiner et all (|2009af ): iGaol (|2009( )). it is clear that such 
events are very rare. So we think this kind of model is 
less likely. 

For GRB 090510, the Fermi collaboration has not 
publis hed the high energy afterglow data yet. Accord- 
ing to lGiuliani et al.l (|2009f ). the high energy photon flux 
recorded by ylG/ii; is - 0.01(V2s)-i-3 cm'^ 8"^ In 

the IC scattering model, the number of the seed photons 
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is needed to be Nsc ^ 10^^. Even all seed photons are in 
near infrared band ('^ 1 eV), the total energy should be 
10^'^ erg, seeming unreasonably large for GRB 090510. 

4.2. Synchrotron radiation of forward shock electrons? 

The spectrum of the synchrotron radiation of 
shocked electr ons can extend to an energy ^ 30 AT/ (1 + 
z) MeV fe.g- ICheng fc WeilfTOOl . where T is the bulk 
Lorentz factor of the emitting region and ~ (1, 2n), 
depending on the comoving acceleration timescale of the 
particles. But usually the IC scattering plays a more im- 
portant role in producing high energy afterglow emission. 
The situation changed in GRB 080319B, the naked-eye 
burst with abundant optical and X-ray afterglow data. 
With the well constrained parameters, Zou et al. (2009, 
Figure 3 therein) have shown that the forward shock syn- 
chrotron radiation dominates over the synchrotron self- 
Compton radiation up to an energy ^ 10 GeV. The de- 
tection prospect for LAT is pretty good. Our estimated 
forward shock parameters of GRB 0809 16C are similar 
to those of GRB 080319B, a strong forward shock syn- 
chrotron GeV emission is nat urally expected (see also 
iKumar fc Barniol Duranll2009[ ). 

In the synchrotron radiation model, the random 
Lorentz factor of electrons emitting > 10 MeV after- 
glow p hotons is so high that rj^^ 1 (^e.g.. lFan fc PiraiJ 
l2006aD . one should take F ~ in calculating Vc oth- 
erwise the radiation flux will be underestimated. For 
GRB 080916C, at t - 400 s i^c < 100 MeV, the flux 
i^ioo MeV = F,^^,^{iyJi^^)-^P~^^/^lOOMeY/hiy,)-P/^ ~ 

9 7 V in-8 Tv plOS /-115n-2 

^ -'^k,55.6*^B,-2.6'^o,-l'^2.6 ^h,29.1 

and the corresponding energy flux is ~ 6.5 x 
10^^ erg cm~^ s""'^, matching the observation 
- 1.2 X 10"^ erg cm^^ g-i p^j. Qj^g 090510, at 
i ~ 5 s, hi/c ^ 18 MeV, so the high energy flux 

i^lOO MeV = F,,„ax(i^c/i^m)-(^-'^/'(100MeV//iZ/e)-P/' 

2.0x10-6 Jy e,^^o.iEi,,54to.\Dl 

.28.26- The corresponding 
energy flux is ~ 5.0 x 10"^ erg cm"^ s~^. The GeV pho- 
ton flux recorded by AGILE is ^ 4 x 10"^ ph cm"^ 
for t - 5 s (see Figure 3 of iGiuliani et all l2009n . 
suggesting an energy flux 6.4 x 10-^erg cm"^ s~^. So 
the observation may be accounted for. 

As shown in iZou et all (|2009[ ) , the synchrotron self- 
Compton radiation of such energetic forward shock peaks 
at TeV energies and may be idea target for the ground- 
based Cherenkov telescopes, like MAGIC (Major Atmo- 
spheric Gamma-ray Imaging Cherenkov Telescope) and 
H.E.S.S. (The High Energy Stereoscopic System). 

5. CONCLUSION AND DISCUSSION 

In this work we have interpreted the high en- 
ergy emission and the afterglow of GRBs G8G916C and 
090510. For the prompt high energy emission of GRB 
080916C with a featureless Band spectrum, the stan- 
dard/unmagnetizcd internal shock model is disfavored. 

^ For the two long bursts E-y > lO^** erg, while for the short burst 
GRB 090510 Ej > lO^^ gj,g ^.11 are at least 1 order of magnitude 



The main reason is that in such a model, the fast shells 
move with a very high bulk Lorentz factor (^ 5 x 10'^) 
and the thermal radiation from their photospheres will 
be too strong to be hidden by the non-thermal emission 
of the internal shocks. As for the idea that the prompt 
GeV photons are the synchrotron radiation of the for- 
ward shock electrons, we predict very unusual X-ray (for 
t < 10"^ s) and optical (for t < 1 day) afterglow light 
curves. The lack of early afterglow observation, however, 
hampers us to test the model. If the prompt GeV pho- 
tons and the soft gamma-rays are from t he same region, 
a non-baryonic c omponent seems needed ( Zh ang fc Pe'eii 
[2009HFa^[200g) . For GRB 090510, the prompt spectrum 
consists of two distinct components. The MeV emission 
may be from the photosphere while the GeV emission is 
produced in the IC scattering of the photosphere photons 
by the shocked electrons. We suggest that the outflow of 
GRB 090510 is baryonic. 

The circum-burst medium of GRBs 080916C and 
090510 is wind-like and ISM, respectively. The stan- 
dard external shock model can reproduce the afterglow 
data reasonably well. The common features are the low 
density of the medium they are expanding into and the 
very high isotropic-equivalent kinetic energy of the out- 
flows. We have proposed a simple method to estimate 
the total number of the seed photons supposing the GeV 
afterglow emission is due to the IC radiation of the for- 
ward shock electrons. Such a model is disfavored be- 
cause the seed photons needed in the modeling are too 
many to be realistic. Though other possibilities, for ex- 
ample the GeV afterglow photons are the synchrotron 
self-Compton radiation of the extended X-ray emission, 
cannot be ruled out, the high-energy afterglow detected 
in these two bursts may be just the synchrotron radiation 
of the forward shock electrons. Our analysis is then in 
support of the "prediction" of Z ou et al. (2009) in GRB 
080319B and the suggestion of Kumar fc Barniol DuranI 
(2009) for GRB 080916C. GRBs 080319B, 080916C and 
090510 are very unusual. They are extremely bright'^ 
and may have a very large initial bulk Lorentz factor. 
Both facts are helpful to give rise to a strong GeV syn- 
chrotron radiation of the forward shock. The number 
density of the circum-burst medium is very low, which 
lowers the detection prospect of the IC radiation compo- 
nent for LAT. For normal GRBs, the detection prospect 
of the GeV synchrotron radiation of the forward shock 
will be much less promising. 
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